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ABSTRACT 

Double-peaked broad emission lines in Active Galactic Nuclei (AGNs) may indicate the existence of a bound 
supermassive black hole (SMBH) binary where two distinct broad line regions (BLRs) contribute together to 
the line profile. An alternative interpretation is a disk emitter origin for the double-peaked line profile. Using 
simple BLR models, we calculate the expected broad line profile for a SMBH binary at different separations. 
Under reasonable assumptions that both BLRs are illuminated by the two active SMBHs and that the ionizing 
flux at the BLR location is roughly constant, we confirm the emergence of double-peaked features and radial 
velocity drifts of the two peaks due to the binary orbital motion. However, such a clear double-peaked feature 
only arises in a particular stage of the binary evolution when the two BHs are close enough such that the 
line-of-sight orbital velocity difference is larger than the FWHM of the individual broad components, while 
the two BLRs are still mostly distinct. Prior to this stage, the velocity splitting due to the orbit motion of the 
binary is too small to separate the emission from the two BLRs, leading to asymmetric broad line profiles in 
general. When the two BHs are even closer such that the two BLRs can no longer be distinct, the line profile 
becomes more complex and the splitting of the peaks does not correspond to the orbital motion of the binary. 
In this regime there are no coherent radial velocity drifts in the peaks with time. Asymmetric line profiles 
are probably a far more common signature of binary SMBHs than are double-peaked profiles. We discuss 
the temporal variations of the broad line profile for binary SMBHs and highlight the different behaviors of 
reverberation mapping in the binary and disk emitter cases, which may serve as a feasible tool to disentangle 
these two scenarios. 

Subject headings: black hole physics - galaxies: active - quasars: general - surveys 



1. INTRODUCTION 

Binary supermassive black holes (SMBHs) a re generic out- 
come s of hierarchical galaxy mergers (e.g., IColpi & Dottil 
l2009l and references therein). A couple of galactic-scale bi- 
nary AGNs were reporte d based on spatially resolved imaging 
and spectroscopy (e.g. , iKomossa et al] 120031 : iBianchi et alj 
l2008HComerford et al.1 l2009bh . These binaries are at pro- 
jected separations of the order of kpc, below which it is dif- 
ficult to spatially resolve both SMBHs at cosmological dis- 
tances. While in principle radio interferometry can resolve 
close binaries down to millarcsecond resolution, it requires 
both black holes (BHs) to be radio sources and so far there 
is only one sub-kpc binary candidate detected in the radio 
(iRodriguez et alJl2006) . 

Characteristic velocity offsets or double-peaked features 
in AGN emission line (both broad and narrow lines) pro- 
files have been invoked to indicate the presence of a binary 
SMBH, even if its spa tial extent is not resol v ed. While this 
idea is not new (ej; ., iHeckman etaP 1 198 It iGaskelll 119831: 
iPeterson et alj Il987l) . it only became feasible recently to 
search for such binary candidates in a systematic way using 
large statistical samples with good spectral quality (most no- 
tably the Sloan Digital Sky Survey' (SDSS) samples). In- 
creasingly larger spectroscopic samples are starting to offer 
a unique opportunity to search for candidate binary SMBHs 
based on t heir spectral properties (e.g., Zhou et al 2004; 
Bonning et al. 2007; Komossaetal. 2008; Comerford et aL 
2p09a; Boroson & Lauer 200 9i: iLiu et anKoiOal: iSmith et ai] 
2009 ; Shields et al. .2009: Wang et alJl2009l: IXu & Komossal 
20091) . This is an important first step towards quantifying the 



viding constraints on hierarchical galaxy merger models and 
predictions for future low-frequency gravitational wave ex- 
periments such as the Laser Interferometer Space Antenna^ 
(LISA). 

While it has become routine to select candidate binaries 
from large spectroscopic data sets, these candidates are less 
secure than spatially resolved cases. The usual difficulty in- 
volves the poorly understood emission line region geometry 
and kinematics even for single black holes (BHs), which may 
mimic a binary system. It is rather difficult to rule out one 
or the other based on a single-epoch spectrum alone. In the 
case of kpc separation binary SMBHs, high spatial-resolution 
imaging and spectroscopy may potentially resolve b oth active 
BHs and therefore confirm the binary nature (e.g, iLiu et alj 
HoiOb). For gravitationally bound binary SMBHs (< 10 pc), 
spectral monitoring and reverberation mapping ma y disen- 
tangle the b i nary scenario from its alternatives (e.g.. | Gaskel]| 
IT9831 lIoTot lEracleous et aT 1997; 'Gezari et all |2007|) . At 
even smaller separations, Loeb (2010) suggested that the time 
dependence of the broad lines due to orbital motion can be 
detected in binaries on the verge of entering the gravitational 
wave dominated inspiral. 

Given the importance of binary SMBHs for galaxy forma- 
tion models and for future low-frequency gravitational wave 
detection experiments, it is crucial to understand the dynam- 
ics and geometry of emission line regions in a binary sys- 
tem and to identify such binaries in a sys tematic way. How- 
ever, despite of the ongoing efforts (e.g.. 'Escala et al.' l2005l: 
Dotti et al. 2006, 2007; Mayer et al. 2007; Bogdanovi d et alj 
20081: ICuadra et al.ll2009l) . it is still challenging to explore the 



frequency of binary SMBHs at various separations, and pro- parameter space in detail in numerical simulations with the 



' http://www.sdss.org/ 



^ http://lisa.nasa.gov/ 
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adequate dynamical range and the necessary input physics. 
On the other hand, observational constraints on the dynamics 
and geometry of AGN emission line regions have not yet con- 
verged to provide reliable inputs for numerical simulations. 
Here we use simple toy models to study the feasibility of us- 
ing broad line diagnostics to identify bound binary SMBHs, 
taking into account both dynamical and ionization effects of 
the two BHs on the combined BLRs. In ^2.1l we develop toy 
models to predict broad line profiles under various circum- 
stances. In ^we discuss the temporal properties of the AGN 
spectrum in different scenarios. Finally, we discuss our re- 
sults in ^ Unless otherwise stated, we will use H/3 as the 
fiducial emission line since this is the best-studied broad line 
in reverberation mapping studies, from which characteristic 
BLR properties are best determined. 

2. BROAD EMISSION LINE PROFILES 

We consider the broad line profiles for bound binary 
SMBHs where the dynamics is dominated by the gravitational 
potential of the two BHs, assuming both BHs are active and 
have their own BLRs. We also assume that the two BHs are 
corotating on a circular orbit. We are interested in binary 
SMBHs with comparable masses (0.3 < C = M1/M2 < 1), 
since binaries with extreme mass ratios are either difficult to 
detect with broad line diagnosis (if line emission is propor- 
tional to BH mass), or difficult to form due to the extended 
dynamical friction time of the galaxy merger. For simplic- 
ity, we also assu me a fixed Eddington ratio AEdd = 0.1 (e.g., 
IShen et ani2008l) . but we note that in practice the two active 
black holes could have different Eddington ratios. For demon- 
stration purposes we will only show binary examples with an 
edge-on view (Figs.[T]|9l). The effect of inclination / is to re- 
duce the radial velocities by sin/. We use to denote the 
binary orbital phase with = 0° when the radial velocities of 
the two BHs are maximal. 

If both BHs have their own distinct BLR, and each BLR 
is corotating along with its BH in the binary orbit, we ex- 
pect to see two time-varying velocity components in the broad 
line profile. The velocity splitting of the two components de- 
pends on the binary orbit separation and the BLR sizes. When 
the orbital separation is large, the velocity splitting of the two 
broad line components is small compared with the broad line 
width and the two components will blend with each other in 
the spectrum. When the two BHs come closer, the velocity 
splitting gets wider, and a double-peak profile may emerge. 
When the two BHs come even closer, both BLRs are dynam- 
ically influenced by the two BHs, and so some BLR material 
becomes circumbinary, leading to a more complex velocity 
structure in the combined line profile. The velocity peaks in 
the broad line, if any, will not simply correspond to the or- 
bital motion of the binary in this case. Eventually, the two 
BHs will get so close that they dynamically affect the BLR 
clouds like a single BH, and the broad line profile may be- 
come single-peaked again. These simple arguments suggest 
that a clear double-peaked broad line profile may only arise 
during a particular stage of the binary orbit evolution, where 
the orbital velocity of the binary is large enough to split the 
two peaks and both BLRs are still mostly distinct. 

Denoting the separation of the binary as d, and the 
BLR radius in a single BH as Rblr, the criterion for both 
BLRs to be bound to their own BHs is Rblr 0.5d. For 
a single BH, assuming the BLR is photoionized by the 
continuum luminosity from the central BH, the Rblr - L 
relation (assuming Rblr oc L'/^) derived from reverbera- 
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Fig. 1 . — Top: Distributions of BLR clouds (projected onto the binary or- 
bital plane) for a binary of two 10** Mq BHs with a separation d = 0.1 pc. 
For clarity we only show a small fraction of randomly selected test particles. 
Different colors indicate BLR clouds initially associated with the two BHs. 
The rotation of the binary is counterclockwise. The observer is located in the 
xy plane at )• = +00, and the radial velocities of the two BHs are maximal at 
this phase. The few clouds that become mixed with the other BH were ini- 
tially on highly eccentric orbits with large apocenter or on hyperbolic orbits 
which travel to the vicinity of the other BH later. We did not remove such 
clouds in our simulation as they have essentially no effect on the derived line 
profile. Middle: Line profile when the radial velocities of the two BHs are 
maximal (orbital phase angle </< = 0°). The dotted lines are individual con- 
tributions from the two BHs. Bottom: Line profile after 1/8 of the orbital 
period (</> = 45°). 
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Fig. 2. — Distributions of BLR clouds (projected onto tlie binary orbital 
plane) and line profiles for a binary of two 10* Mq BHs with a separation 
d = 0.05 pc. Notations are the same as in Fig.[T] The cyan dashed line in the 
middle panel shows a disk emitter model (see the text for details). 
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Fig. 3. — Distributions of BLR clouds (projected onto the binary orbital 
plane) and line profiles for a binary of two 10* Mq BHs with a separation 
d = 0.02 pc. Notations are the same as in Fig.[T] 
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tion mapping (e.g.. iKaspi et alj l2000t iPeterson et alj l2004t 
iMcLure & Jarvisll2002HKaspi et al.ll2007h~ 



gives 



: 2.2 X 10-2 r^') 



1/2 



1/2 



pc , (1) 



where M, is the BH mass and Andd is the Eddington ratio of 
the broad Hne AGN^*. This T^BLR-i- scaHng impHes that the 
ionization parameter, U oc L/r^, is roughly constant for a par- 
ticular line species (assuming constant electron density). The 
orbital period for a circular orbit at this location is: 



'orb 



= 30 



0.1 



1/4 



yr. 



(2) 



Assuming the BLR is virialized, the FWHM of a single 
broad component is: 



VpWHM = 



/ GM, 



KfRB 



1/2 



^4200 



^Edd 

"oT 



-1/4 



VlO^Mc 



1/4 



J_ 

1.4 



-1/2 



kms 



(3) 



where / is the virial coefficient accounting for our ignorance 
of the BLR geometry, and is of order of unity. A recent deter- 
mination of / based on the reverberation BH masses a nd BH 
masses from the M, - a relation gives / « 1 .4 ( Onke iTet al.l 
12004)'*. 

If the two BLRs are distinct and simply corotate along with 
their own BHs in a binary system, the line-of-sight (LOS) ve- 
locity splitting of the two broad line peaks changes periodi- 
cally in time f : 



Vlos = 



GM,, 



1/2 



:6300 



M, 



sin/sin(27rf/P) , 

1/2 



X sin/sin(27rf/P) kms"' 



-1/2 



(4) 



where Mjot = Mi +M2 is the total mass of the two BHs, / is 
the inclination of the binary orbital plane relative to the LOS, 
and 



P = 27rc/^/2(GMtot)-'/2 

ipcy V 10*^^0 



= 9.5 X 10- 



M„ 



-1/2 



yr 



(5) 



is the orbital period of the binary. 

From these simple calculations ([B-® we can investigate 
how well-separated the double components could be without 
violating the assumption that each BLR is mostly under the in- 
fluence of only one BH. In order to see a clear double-peaked 

^ The RblR"^ relation is usually calibrated using the restframe 5100 A 
continuum luminosity. We have adopted a bolometric correction 1^^^^^^^^ = 

10 to convert the continuum luminosity to bolometric luminosity. 

Note that we have used a different definition of the virial co efficient from 
others . The coefficient / here con'esponds to the coefficient e in lOnken et alj 
<2004h . 



feature we require vios.max ^ i'fwhm, while in order to have dis- 
tinct BLRs we need /?blr ^ d/2. Assuming ^ = M1/M2 < 1 
these two criteria imply: 

1/2 /I / fl^ \ 1/2 



0.044 



M2 



lO^M, 



pc<d< 0.063 



1-1- A/ M2 



V lO** Mr: 



pc . 

(6) 



Therefore, there is only a narrow window of separation (as 
well as mass ratio ^) within which a double-peaked broad line 
profile may emerge without violating the assumption that both 
BLRs are distinct. 

2.1. Model Setup 

Next, we investigate the broad line profile for a binary 
SMBH in a more quantitative manner. First we describe our 
approach to model the BLRs of a binary SMBH. 

We start with a simple prescription for the BLR around a 
single BH, where the BLR is assum ed to be an assembly of 
discrete clouds (e.g.. iPetersonI [19971) . Then two BHs with 
their individual BLRs are placed on a circular orbit. We in- 
tegrate numerically the orbits of individual clouds (treated as 
test particles) in the circular restricted three body problem. By 
imposing an ionization condition (see below) we identify the 
clouds that will radiate the line emission and determine the 
line profile based on the LOS velocity distribution of these 
line-emitting clouds. This procedure is detailed below. 

First, we need to specify a model for the BLR around a 
single BH. Despite decades of research, the detailed structure 
of BLR is still poorly constrained. The most powerful ob- 
servational tool to study BLR structure is reverberation map- 
ping. But even with the best-studied reverberation mapping 
sample, there is still no complete conse nsus on the general 
BLR structure (e.g.. iDennev et al]|2009l) . Nevertheless, re- 
verberation mapping does provide a characteristic scale for 
the BLR radius, given by the Rblr-L relation in Eqn. ([T]). 
In a handful of cases, reverberation mapping of different line 
species in the same system shows that lines with larger ioniza- 
tion parameters have sma ller radii from the central BH (e.g., 
iPeterson & Wandeill2000h . and that the velocity inferred from 
the line width is consistent with virialized motion. 

Motivated by these observations, we adopt the following 
simple model for the BLR around a single BH. For a sin- 
gle BH, we uniformly populate clouds as test particles within 
a spherical shell with inner and outer radii VoTS and \/L2 
of the radius in Eq. ([T]), i.e., the flux or ionization parameter 
(assuming constant electron density) required to photoionize 
a particular line species in the BLR clouds is roughly con- 
stant within ±20% (the exact value of the percentage does not 
change our conclusions). The velocities of those clouds are 
generated from a Maxwellian distribution whose ID disper- 
sion is determined from the virial relation, a = ^ GM,/{3r), 
with random orientations. Since it is unphysical to restrict 
these clouds within a perfectly thin shell, we evolve the Ke- 
plerian orbits of these clouds in the single BH system using 
standard analytical formulae (e.g. Murrav & Dermott 1999). 
The system quickly establishes a quasi-steady state configu- 
ration where the spatial and velocity distributions converge 
after ^ foib (Eqn. |2]). The resulting quasi-steady-state distri- 
butions are taken as the initial cloud distributions around a 
single BH^. These clouds have a radial distribution peaked 
around the initial shell location but also extend to larger and 

' We do not remove the small fraction (~ 10%) of clouds on hyperbolic 
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smaller radii, and the Gaussianity of their velocity distribution 
is mostly preserved. We generate 5000 test particles for each 
BH. As discussed later, the sporadic snapshots of these test 
particle trajectories will be combined to compute the cloud 
distributions in a binary system. The resulting line profiles 
are very smooth even though the total number of test particles 
used is modest. 

This simple model for the BLR around a single BH is highly 
idealized, given the limited understandings of actual BLR 
properties. It nevertheless reproduces line profiles and char- 
acteristic BLR sizes that are consistent with observations. We 
have assumed that only clouds at a narrow range of distances 
from the BH can be line-emitting and constitute the BLR. 
In reality the BLR in some AGNs could have a larger spa- 
tial extent, as inferred in some (albeit limited) reverb eration 
mapping studies (e.g., Dennev et al. 2009; Horne et al 120041 
and references therein). This can be understood if the elec- 
tron density in the BLR varies with radius, so that the flux 
(oc r"-) required to produce the proper ionization parameter 
also has a larger range in photoionization models. However, 
given the smal l scatter in the mean R^lr-L relation (e.g., 
iPetersonllMoh . our fiducial choice of a relatively thin shell 
geometry for the BLR should be a reasonable approximation 
at least for the majority of the line-emitting clouds in the BLR. 
One caveat is that we are not including a significant fraction 
of clouds that are either far more distant or closer in than the 
characteristic BLR distance. Such clouds will not be line- 
emitting in the single BH system because the ionization pa- 
rameter is either too high or too low. Clouds much closer in 
will be tightly bound to their BH and will not contribute to the 
line emission in the binary system either. However, if there 
are a large amount of cold clouds orbiting outside the BLR in 
single BHs, these clouds are easier to become circumbinary 
and may be ionized by the combined continua from both BHs 
in the binary system, thus making the double-peaked feature 
less prominent. Below we proceed with our fiducial model for 
the BLR around single BHs, and the effects of such a distant 
cloud reservoir will be further discussed in Sec 14. II 

We combine the two BH plus clouds systems to form a 
binary on a circular orbit. In the frame corotating with the 
binary, we derive the instantaneous locations and velocities 
of these clouds by orbital integrations of the restricted three- 
body problem. We are only interested in the temporal be- 
haviors over a period that can be monitored on human life 
timescales, so we integrate the system for a few hundred years 
sampled with 1,000 evenly distributed temporal snapshots. 
Long-term stability of the cloud orbits are not consid ered in 
this p aper^. We use a Bulirsch-Stoer integrator (Pre ss et al.l 
1 19921) and the relative accuracy tracked by the Jacobi constant 
was well below 10"'" over the course of the integration. We 
then combine all the snapshots as a quasi steady-state config- 
uration of clouds in the binary system, so as to improve the 
statistics of the cloud distributions. We found that the cloud 
distributions (in real and velocity space) in individual snap- 
orbits in the initial distributions, since such clouds will either be bound in 
the binary system, or will escape to large distances and not contribute to line 
emission anyway. We veritied that removing these clouds does not have any 
effect on the derived line profiles. 

^ It is possible that some of these orbits are unstable over longer timescales 
and the long-lived orbits may have somewhat different spatial and velocity 
distributions. On the other hand, the individual BLRs are likely non-static 
and some mechanisms may exist to replenish the clouds continuously (such 
as a wind from the accretion disk). It is beyond the scope of this paper to take 
into account these complications, but we note that our integration time may 
be too short for some of the cases studied here. 



shots are similar (but with much poorer statistics) to those 
from the combined snapshots, after the initial ^ 10% of in- 
tegration time, i.e., the velocity distribution of BLR clouds 
quickly establishes a quasi-equilibrium in much less than one 
binary orbit time. Although clouds are still constantly ex- 
pelled by the binary afterwards, the loss of clouds does not 
have significant effects on the velocity structure of the BLR 
over the course of the integration. 

Once we have the locations and velocities of the clouds, 
we determine line-emitting BLR clouds under the assumption 
that clouds are illuminated by both BHs, and that the com- 
bined flux (or ionization parameter) required to ionize a BLR 
cloud is roughly constant (within ±20%, assuming a constant 
electron density). Quantitatively, the locations of these line- 
emitting clouds satisfy the following constraint: 

/AEdd.Iw Ml s /AEdd,2w M2 ^ 

0.8 ^ ^ 0.1 ^''10'^Mq^ 0.1 'HO^Mq' 
(2.2 X 10-2)2 - (ri/pc)2 (r2/pc)2 



- (2.2 X 10-2)2 ' ' ' 

where ri and r2 are the distances to the two BHs. This con- 
dition follows the photoionization criterion for single BHs in 
Eq. dB. 

Based on the procedure outlined above, we derive the broad 
line profile for binary SMBHs at any orbital phase and hence 
the temporal variation of line profile due to the orbital motion 
of the binary. 

Figures [l][3] show several examples for an equal-mass 
binary with Mi = M2 = 10** Mq, and separations of li = 
0.1,0.05,0.02 pc, with integration time 2, 6 and 22 times the 
binary orbital period given by Eq. (|5]l, respectively. Only the 
intermediate separation case (d = 0.05 pc) satisfies the con- 
straint of Eq. (|6]l. As the two BHs come closer, their BLRs 
become less distinct, and eventually the two BHs start to af- 
fect both BLRs in terms of their dynamics and illumination. 
In the example with the closest separation (e.g., Fig.O, there 
is practically a single BLR surrounding the two BHs, and the 
velocity structure of these BLR clouds are more complex than 
in the distinct BLR cases. Many BLR clouds are now on 
circumbinary orbits (such as horseshoe or tadpole-type orbits 
around the triangular Lagrangian points) and there are no co- 
herent radial velocity shifts in the two peaks due to the orbital 
motion of the binary. 

Figures |4l|6] show several examples for a non-equal-mass 
binary with Mi = 3 x 10^ Mq and M2 = 10** Mq and separa- 
tions of d = 0.1,0.05,0.02 pc, where the BLR clouds associ- 
ated with each individual BH are given a flux weight propor- 
tional to the BH mass. The integration times are 2, 6 and 22 
times the binary orbital period, respectively. Again, only the 
intermediate separation case {d = 0.05 pc) satisfies the con- 
straint of Eq. (|6]l. The qualitative difference from the equal 
binary case is the increase in line asym metry and decrease in 
the prominence of the double peaks (see iPopovic et al.ll2000l 
for some similar line profiles generated with different model 
prescriptions for the binary BLRs). 

2.2. Disk Emitters 

An alternative interpretation of double-peaked broad line 
profiles is the disk emitter scenario, where the anoma- 
lous broad line emission originates f rom a rela t ivistic 
accretion disk around a single BH (IChen et al.l 119891; 
lEracleous & HalpernI IT994t lEracleousetal.lll995l) . In fliis 
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Fig. 4. — Distributions of BLR clouds (projected onto tlie binary orbital 
plane) and line profiles for a SMBH binary with Mi = 3 X lO' Mq and M2 = 
10** Mq with a separation J = 0.1 pc. Notations are the same as in Fig.[T] 

case the blueshifted and redshifted components originate from 
the part of the disk moving towards and away from the 
observer. The disk emitter model has been successful in 
reproducing the line profi le in many double-peaked broad 
Une AGNs (e.g IChen et al.a989;,Eracleous & Halperal994t 
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Fig. 5. — Distributions of BLR clouds (projected onto the binary orbital 
plane) and line profiles for a SMBH binary with = 3 X lO' Mq and Mi = 
10** Mq with a separation d = 0.05 pc. Notations are the same as in Fig.[T] 
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Fig. 6. — Distributions of BLR clouds (projected onto tlie binary orbital 
plane) and line profiles for a SMBH binary with Mi = 3 X lO' Mq and Mi = 
10** Mq with a separation d = 0.02 pc. Notations are the same as in Fig.[T] 



|2009|) . As an example, in the middle panel of Fig. |2] we show 
a disk emitter model for an inclined (inclination I = 25°) el- 
liptical disk (with eccentricity e = 0.1) with inner and outer 
radii of 500rg and SOOOri, (where = GM/c^ is the gravita- 
tional radius), specific intensity c>c r"^, an internal turbulent 
broadening a = 800 kms ~' and a major axis or ientation of the 
elliptical disk 4>q = 60° (lEracleous et alJll995h . around a BH 
with M = 10^ Mq. The disk emitter model clearly shows a 
double-peaked profile, resembling that for a SMBH binary. 
But the temporal variations of the broad line in the disk emit- 
ter scenario are different from those for the binary SMBH sce- 
nario, as we discuss next. 

3. TEMPORAL VARIATIONS 

3.L Radial Velocity Drifts 

A definitive signature of a binary SMBH is the time drift in 
the radial velocities of the decomposed two components, as 
resulting from the orbital motion of the two BHs. Unfortu- 
nately, the typical orbital time is usually much longer than a 
few years, and in order to detect radial velocity drifts of the 
double peaks the two BLRs must be distinct. From the previ- 
ous sections we know that the optimal configuration to detect 
such a binary is when the two BLRs are just touching each 
other, such that the two BLRs are still mostly distinct while 
at the same time the velocity splitting of the two components 
is larger than the line width. In the case of two equal mass 
BHs, substituting d = 27?blr in Eq. (|5]l where 7?blr is given 
by Eq. Q, we have fopt 62(M,/10'* Mq^I^ yr This means 
it is less challenging to detect radial velocity changes in low 
mass SMBH binaries by spectral monitoring. In the case of 
two M, = 10^ Mq, the time span between double-peaked and 
single-peaked profile is only /'opt/4 ^ 5 yr. The limitation of 
low-mass binaries is that they cannot be easily observed out 
to high redshifts due to their relatively low luminosities. Nev- 
ertheless, these low-mass systems provide good test cases in 
the nearby Universe. 

Figure [7] shows changes in the overall line profile due to 
the orbital motion of the binary for an intermediate separa- 
tion case (the example in Fig. |2]) and for a close separation 
case (the example in Fig. |6]l. While the radial velocity drifts 
of the two peaks in the intermediate separation case are ap- 
parent, there are no coherent drifts in the radial velocities 
of the two peaks in the close separation case when the two 
BLRs are no longer distinct (see the upper panel of Fig. |6l). 
The non-detection of coherent radial velocity drifts in some 
of the double-peaked broad line AGNs may t hen rule out the 
existence of two d istinct BLRs (e.g., Eracleous et alj IT997I : 
iGezari et alJ l2007h . but cannot rule out the possibiUty of a 
close SMBH binary surrounded by a circumbinary BLR, as 
pointed out in lEracleous et alj (Il997h . Therefore additional 
tests are required to distinguish the binary and disk emitter 
scenarios in these cases^ (see below). 

3.2. Reverberation Mapping 

A better way to distinguish the binary scenario and 
the disk emitter scenario is reverberation mapping^ (e.g., 

^ In addition to the reverberation mapping method discussed in ^3.21 there 
could be variations due to the orbital motion (Doppler effect) of the contin- 
uum emitting region of the accretion disk around each BH (Kocsis & Loeb 
2010, in preparation). 

^ In principle, one could argue that independent variations in the two veloc- 
ity components of the broad line may be sufficient to distinguish the binary 
interpretation from disk emitters (e.g., IGaskelll2010h . However, reverbera- 
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M^ = M2=10° M3^,, d = 0.05 pc Mi = 3x1 0'M3^,,M2= 1 0°M3^,,d = 0.02pc 




V|os (kms ^) V,„3 (kms"^) 

Fig. 7. — Time series of the radial velocity drifts in the line profile for an edge-on binary SMBH and for 1/4 of the binary orbital period. Left: a binary with 
Ml = M2 = 10^ MQ,d = 0.05 pc, and P 75 yr (see Fig.|2). The two BLRs are distinct in this case, and the drifts in radial velocities of the two peaks are 
apparent. Right: a binary with Mi = 3 X lO' Mq, M2 = 10* Mq, d = 0.02 pc and P 23.5 yr (see Fig.|6). The two BLRs are no longer distinct in this case, and 
there are no coherent drifts in the radial velocities of the two peaks. 
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V|os (kms"') V,„3 (kms"') 

Fig. 8. — Responses of broad line profiles to continuum variations of one BH for the wide separation case. Shown here is an example of two 10* Mq BHs with 
a separation d = 0.1 pc, emitting at Andd =0.1 and in an edge-on view, with orbital phases and 7r/4. Solid lines show the overall line profile and dotted lines 
show individual components. The black lines show the original line profile, and the colored lines show the final profile after the reverberation completes. The 
bottom panels show the velocity-resolved flux changes. Even though the overall line profile is only mildly asymmetric, the fact that only one BLR is reverberating 
is seen once we decompose the broad line profile into two Gaussian components. 
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Fig. 9. — Left: Responses of broad line profiles to continuum variations of one BH for a binary with an intermediate separation. Shown here is an example 
of two 10* Mq BHs with a separation d = 0.05 pc, emitting at \Edd =0.1 and in an edge-on view, with orbital phase = 0. We assume a 40% increase in the 
continuum of one BH and derive the line response at different late times. Solid fines show the overall line profile and dotted lines show individual components. 
The black lines show the original line profile. The fine response is completed at the last shown time epoch (r = 90 days). The continuum variability of one 
BH also has effects on the BLR clouds of the other BH, which are delayed compared to the response of its own BLR clouds and are less prominent. Right: 
Reverberation mapping for a disk emitter around a single 10* M^BH. The amplitude of the continuum increase is 40%. The black line shows the original line 
profile. The disk model is the same as in the middle panel of Fig.12] The line response is completed at the last shown time epoch (r = 30 days). 
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iPeterson etallligS?!: iGaskellllloToh . Due to the different ge- 
ometries of the BLR in the two scenarios, the response of the 
line to the variations in the continuum follows different pat- 
terns. 

In the case of wide separation binaries, each BLR is illumi- 
nated by its own BH, and therefore only responds to the lu- 
minosity variations of its BH. We demonstrate this case with 
the example of two 10^ Mq BHs with a separation 0.1 pc. 
The orbital time of the binary is ^ 200 yr, while the typi- 
cal variability timescale for these BH masses is a few months 
(e.g., P eterso n 1997). Once the continuum of one of the BHs 
varies, it takes ~ 1 month for the associated BLR emission 
to vary, and the orbital motion of the binary during this light 
travel time is negligible. The effects of continuum variations 
of one BH on the overall broad line profile are shown in Fig. [8] 
for several variability levels and two orbital phases, after the 
broad line has completed its reverberation {t > 35 days). Since 
the binary separation is large, the broad line is not double- 
peaked. But the line response introduces an asymmetry that 
correlates with the amplitude of the continuum variation. This 
asymmetry is more prominent at smaller orbital phase angles 
where the two velocity components overlap less in the spec- 
trum. 

In the case of intermediate separation binaries, the two BHs 
are close enough such that part of the BLRs are illuminated 
by both BHs. Consider the example of two 10** Mq BHs 
with a separation 0.05 pc, i.e., about twice the size of a single 
BLR. In this case, we use the configuration of BLR clouds 
derived in our orbital integrations to compute the line pro- 
file changes due to continuum changes of one BH. The or- 
bital time of the binary is ^ 75 yr, hence the orbital motion 
is negligible compared to the light travel time to the BLR, 
which is a few months. However, due to the proximity of the 
two BHs, the continuum variation in one of the BHs also af- 
fects the other BLR, although the effect is delayed and less 
prominent than those on its own BLR. The left panel in Fig- 
ure |9] demonstrates these effects for an edge-on view and or- 
bital phase = 0°. For comparison, the right panel shows the 
expected line responses to continuum variations for the same 
disk emitter model shown in Fig.[2]at various times, where the 
line response is complete at the last epoch. Although the high- 
velocity wing responses more quickly than the central part of 
the line (because the former emission originates closer to the 
BH), the two peaks change simultaneously. Even though ini- 
tially the broad line profiles are similar in the binary case and 
in the disk emitter case, their subsequent line responses are 
different: ii) the disk emitter completes the response of the 
line much faster than the binary due to the proximity of the 
disk to the central BH; ii) both peaks are varying in propor- 
tion to each other in the disk emitter case, while in the binary 
case one peak has larger amplitude changes than the other; 
and iii) the two peaks reverberate simultaneously in the disk 
emitter case, whereas in the binary case one peak reverberates 
faster than the other 

In the case of close separation binaries where the two BLRs 
are already merged, the behavior of line response to contin- 
uum variations of one of the BHs is similar to that of the in- 
termediate separation binary, but the relative increase in the 
two peaks is less discrepant than in the intermediate separa- 
tion binary. Nevertheless, reverberation mapping will be a 
useful test here because the tests on radial velocity drifts are 
inconclusive, as discussed in ^3.11 

4. DISCUSSION 
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Fig. 10. — Results for an extended cloud distribution around a single M = 
10^ Mq BH. Upper: Radial distributions of clouds around the single BH as 
initial conditions for the numerical orbit integr ation s. The green dashed line 
is the fiducial cloud distribution described in j|2.1l which peaks ai'ound the 
characteristic radius giv en b y Eqn. Q. The black line is a more extended 
cloud distribution (see j|4.1l for details). Bottom: The resulting line profile 
for the extended cloud distribution and for an equal-mass binary (Mot = 2 X 
Mq) with a binary separation d = 0.05 pc at phase angle (/> = 0°. The two 
line components ai'e much less distinct than in the previous case (e.g., Fig. 
|2) because a larger fraction of the outmost BLR clouds are on circumbinary 
orbits. This tends to dilute the distinction between the two components from 
BLR clouds orbiting around individual BHs. 



4.1. Complications in the Realistic Situation 

Our simple prescription for the BLR of a SMBH binary is 
by no means realistic, especially for the closest separation 
cases studied here, where there are no longer two distinct 
BLRs. The BLR models we adopted are close to models in 
which the orbits are random, and the dynamics is dominated 
by the gravitational potential from both BHs. The novelty of 
our approach is to include the effects of the two BHs in terms 
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Fig. 1 1. — An example of double-peaked broad line AGNs from SDSS. 
The black lines show the original SDSS spectrum (upper) and the contin- 
uum/iron flux subtracted spectrum (bottom). The red and magenta lines show 
the overall model fits. The orange and blue lines ai'e the power-law contin- 
uum and iron template fits. The three Gaussian components in cyan are the 
narrow lines H/3, [O III] A4959 and [O III] A5007, fixed to have the same 
redshifts and line widths. The two green Gaussian components are for the 
double-peaked broad line profile, which are separated by ^ 3000 kms"' , and 
have FWHMs ~ 2200 kms"' and ~ 10000 kms"' for the blueshifted and red- 
shifted components respectively. This object and many others will be good 
candidates for spectral monitoring and reverberation mapping programs. 

of both the clouds dynamics and illumination. Our treatment 
is more quantitative than earlier qualitative arguments that the 
double-peake d components are from individual BLRs (e.g., 
lGaskellll2010i and references therein). 

In our BLR model for single BHs we assumed a rather sim- 
plistic thin shell distribution of clouds. To check the sensitiv- 
ity of our results to this assumption, we examine the effects 
of a more extended cloud (not all line-emitting clouds) dis- 
tribution for single BHs. The upper panel of Fig. [TOl shows 
an example for the distribution of clouds around a 10*^ Mq 
BH, where the clouds were initially populated between 0.5 
and 2 times the characteristic BLR size with a power-law 
number density n(r) cx r"' and then relaxed for 30 years us- 
ing their Keplerian orbits. The initial random velocities of 
each cloud are assigned using the same scheme described in 
92.11 The starting configuration for the numerical orbit inte- 
gration is a more extended distribution of clouds compared to 
our fiducial model. We integrate an equal-mass binary sys- 
tem (Mtot = 2 X IO^Mq) using the new single BLR model and 
compute emission line profiles in the same way as in 92.11 
The bottom panel of Figure[TO]shows the resulting line profile 
at phase angle = 0° and with a binary separation li = 0.05 pc. 
Compared with our previous result, i.e., the middle panel of 
Fig. |2] the double-peaked feature is much less prominent. 
This is expected because now the two cloud regions over- 
lap more than in the previous case and more clouds become 
circumbinary, diluting the distinction between the two BLRs. 
In practice, the two emission line components will not have 
equal strength, so some asymmetry is expected in the overall 
line profile, similar to those shown in Fig.|5] Secular changes 
in the line profile due to the orbital motion of the binary, or 
velocity resolved reverberations of the blue and red wings of 
the line, can still be monitored even though the two peaks are 
blended with each other. 

Our toy models confirm the feasibility of using spectral 
monitoring (for radial velocity drifts) and reverberation map- 
ping (for line responses) to disentangle SMBH binaries from 
disk emitters. These observations can also be achieved with 
less expensive narrow-band filters. However, there are some 



practical difficulties with these techniques. The spectral mon- 
itoring for radial velocity changes works best for binaries with 
d « /^BLR.i +^BLR.2- At larger binary separations, the two 
broad components will blend with each other in the spec- 
trum, making it less likely to be flagged as a binary candi- 
date, and the orbital period is too long to be detectable. At 
smaller separations, the BLRs are no longer distinct and the 
velocity structure becomes more complex with no coherent 
radial velocity drifts in the peak locations. Thus, the spec- 
tral monitoring is suitable for identifying low mass SMBH 
binaries (M ^ 10^ Mq). On the other hand, reverberation 
mapping is a powerful tool for distinguishing a SMBH binary 
from a disk emitter, but more complex BLR geometries and 
kinematics (such as those involving infl ows or outflows) will 
certainly cornplicat e the situation (e.g.. ISergeev et alj 119991: 
iDenney et al]l2009h . 

4.2. Case Studies: 3C 390.3, SDSS J 1536+0441 and others 

Although there are a few dedicated spectral mo nitoring pro- 
grams for double-peaked broad line objects (^e.g..' Gezari et alJ 
i2007, and references therein), there is currently no reverber- 
ation mapping program for a large sample of such objects. 
Among the ^ 40 AGNs that are includ ed in the reverberation 
mapping sample (iPeterson et al.ll2004h . there are several ob- 
jects that clearly show double-peaked broad line features. In 
particular 3C 390.3 is a stro ng double-peak object with good 
reverberation mapping data (iDietrich et al.llT9 98). The time- 
ordered data of this source show that the blueshifted and red- 
shifted components respond to the continuum variations al- 
most simultaneously. Thus, it is more likely that the double- 
peaked emission originates from a disk rather than from two 
corotating BLRs in a binary system. 

The quasar SDSS J1536-I-0 441 was recently susp ected 
to be a sub-pc binary SMBH (iBoroson & Laued 12009 ) be- 
cause of its double-peaked broad H/3 line in the SDSS 
spectrum. It has therefore received much attention (e.g., 
Chornock et al. 2010; Gaskell 2010; Wrob el & Laoj |2()(M 
iDecarli et al..,2009: .Tang & Grindlay, ,2009.) . The discovery 
of an additional redshifted component, most nota bly in Ha 
(e.g.. ' Chornock et al.ll20Tol: iLauer & Boroson"2009), favored 
a disk emitter origin for the double peaks rather than a binary 
SMBH. If we nevertheless assume this is a binary SMBH and 
use constant Eddington ratios AEdd = 0.1 for both BHs, and 
the FWHM values measured in Bo roson & Laued (l2009h . we 
get from equation Mi = 8 x 10** Mq and M2 = 2 x 10'' Mq 
for the red and blue systems, respectively, similar to those 
reported in Bor oson & Laued (l2009h . We also derive BLR 
sizes Ri « 0.063 pc and R2 « 0.01 pc. Substituting the 
BH masses and vios = 3500 kms"' in equation (HI we get 
d = 0.27[sin/sin(27rf/P)]^ pc. However, in order to produce 
the comparable strength of both components and hence two 
distinct peaks, it requires that the smaller BH (blue compo- 
nent) is ~ 40 times more efficient at producing the broad- 
line emission than the larger BH. Fur thermore, by compar- 
ing the spectrum taken ^ 1 yr after dChornock et alJl2010l) 
with the original SDSS spectrum, it appears that both the 
blueshifted and the redshifted components become slightly 
weaker whereas the central component (which would be the 
classic BLR in the disk emitter scenario) remains almost un- 
changed. This strengthens the association of a disk emitter 
origin with the double-peaked feature. 

There have been a significant number of double-peaked or 
highly asymmetric broad line AGNs known in the literature, 
and large spectroscopic surveys such as SDSS are provid- 
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ing many more (e.g.. IStrateva et alj 120031: IShen et al.ll2010l) . 
An example, SDSS J0012-1022'', is shown in Fig.[IT] whose 
broad line profile can be well fitted by two components. Note 
that the velocity splitting of this object is smaller than the 
FWHM of the redshifted component, and both broad compo- 
nents are reasonably well-fitted by a single Gaussian. If it is of 
a binary origin, its separation must be large enough such that 
the two BLRs are mostly distinct. Even though its profile is 
not obviously double-peaked, it still stands out as a promising 
binary candidate. It would be interesting to perform spectral 
monitoring and reverberation mapping for a large statistical 
sample of such objects in order to determine their nature. 

4.3. The Frequency of Spectroscopic Binaries 

There is currently some tension between theoretical expec- 
tations and observations on the frequency of SMBH binaries: 
if major mergers of gas-rich galaxies is the triggering mech- 
anism of quasar activity, then the expected binary fraction is 
very high; on the other hand, the observed frequency of bi- 
nary SMBHs is less than a few percent beyond pc scales, and 
much lower on pc to sub-pc scales. We can parameterize the 
"observable" binary fraction at different stages as the product 
of several factors: 

/obs = /bin X /active X /geo X /tech X (Tphase/Tigso) , (8) 

where /bin is the intrinsic binary fraction (/bi„ w 1 in the 
merger hypothesis for quasar activity), /active is the probability 
of both BHs being active, /geo is the observable fraction due 
to orientational effects (inclination and orbital phase), /ech is 
the observable fraction due to the specific technique used (i.e., 
spectroscopic methods or spatially resolved imaging) which 
depends on the data quality of the observing program and 
folds in the complications associated with the emission line 
region geometry, Tphase is the time span during the specific 
stage of binary evolution, and tqso is the lifetime of quasars. 
All these quantities except for tqso are functions of the evo- 
lutionary stage of the binary. 

' Strictly speaking, this is not a double-peaJced profile, but an asymmetric 
profile with an extended red wing. Nevertheless it was included in the double- 
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In this paper, we have focused on the broad Une diagnosis 
technique, which are for parsec to sub-parsec binaries before 
the binary enters the gravitational wave dominated regime. 
Therefore, we expect /active ~ 1 since in gas-rich mergers a 
nuclear gas disk on pc to 100 pc s cales may form and feed 
both BHs (e.g., iMayer et ani2007h . For the effects of ran- 
dom orientations and orbital phases it is reasonable to adopt 
/geo ^ 0.25, since nearly edge-on systems with radial motions 
are most likely to be detected, /ech is difficult to quantify 
without a dedicated program and related Monte-Carlo sim- 
ulations, but it is unlikely that /ech is close to unity. Bear- 
ing in mind the large uncertainties, the time span during 
this stage can be estimated as ~ 10^-10^ yr if gas drag is 
the dominant mechanism that shrinks the binary orbit (e.g., 
lEscala et al.ll200l iDotti et al.l l2006': Ivla ver et al.ll2O07l) . The 
quasar lifetime is not very well constrained and typical values 
are tqso ^ 10''- 10** yr. Taken together, the fraction of observ- 
able parsec to sub-parsec binaries based on broad line diagno- 
sis is less than a few percent. This is still higher than the fre- 
quency of the known parsec to sub-parsec binary candidates, 
but lower than t he frequency of known double-peaked b road 
line AGNs (e.g.. IStrateva et al.ll2003l: iGezari et alj|2007l and 
references therein). Of course, these are very crude estimates, 
and measuring the actual observed frequency will serve an im- 
portant role of testing theoretical models of binary formation 
and evolution (e.g.. iVolonteri et aL,2009.) . 
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